A semiclassical description of the Autocorrelations in Nuclear Masses 



Antonio M. Garcia-Garci'a 
Physics Department, Princeton University, Princeton, New Jersey 08544, USA and 
The Abdus Salam International Centre for Theoretical Physics, P.O.B. 586, 34100 Trieste, Italy 



O ■ 

o ■ 

(N : 

s . 

<: 

00 1 



> 
o 

00 

o 

\o 
o 



o 

=5 



X 



Jorge G. Hirsch and Alejandro Frank 
Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, AP 70-543, 04510 Mexico DF, Mexico 

Nuclear mass autocorrelations are investigated as a function of the number of nucleons. The fluc- 
tuating part of these autocorrelations is modeled by a parameter free model in which the nucleons 
are confined in a rigid sphere. Explicit results are obtained by using periodic orbit theory. Despite 
the simplicity of the model we have found a remarkable quantitative agreement of the mass autocor- 
relations for all nuclei in the nuclear data chart. In order to achieve a similar degree of agreement 
for the nuclear masses themselves it is necessary to consider additional variables such as multipolar 
corrections to the spherical shape and an effective number of nucleons. Our findings suggest that 
higher order effects like nuclear deformations or residual interactions have little relevance in the 
description of the fluctuations of the nuclear autocorrelations. 
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I. INTRODUCTION 

As a consequence of the strong nuclear interaction, the 
nuclear mass M is not just the sum of the individual 
nucleons. The difference between these two quantities is 
an indicator of the stability of a given nucleus, the larger 
the difference the more stable is the nucleus. An accurate 
description of this binding energy as a function of the 
number of neutrons and protons is a recurrent research 
topic in nuclear physics and nuclear astrophysics Q. 

The semi-phenomenological liquid drop model, in 
which the nucleus is described as a very dense, charged 
liquid drop, is the oldest and simplest approach to this 
problem Q. It provides a qualitative description of the 
binding energy though it fails to capture features related 
to the quantum nature of the single particles (neutrons 
and protons) inside the nucleus. This is clearly observed 
in Fig. ^ where we have plotted the difference between 
measured masses Q and Liquid Drop Model (LDM) pre- 
dictions as a function of the proton number N, mass 
number A, neutron number Z, and as an ordered list 
0, 0] • The sharp valleys and round peaks which remain 
after the removal of the smooth LDM mass contribution 
contain information related with shell effects due to the 
quantum motion of the individual nucleons inside the nu- 
cleus, nuclear deformations, and nuclear residual interac- 
tions. One of the main goals of the present paper is to 
further investigate the details of these corrections. 

Most theoretical descriptions of nuclear mass models 
have as a starting point the general expression 



M = M + SM 



(1) 



where M is a smooth function of the number of nucle- 
ons, usually the liquid drop mass formula. By contrast 
SM is a fluctuating function in the number of nucleons 
which accounts for the quantum nature of protons and 
neutrons within the many body problem. There is a va- 
riety of nuclear mass models in the literature, two of the 



most broadly utilized are the finite range droplet model 
(FRDM) [8j , which combines the macroscopic effects with 
microscopic shell and pairing corrections, including ex- 
plicit deformation effects and the Strutinsky procedure 
Q , and, on the other hand, the Dufio and Zuker (DZ) [Io| 
model, where the microscopic corrections are functions of 
the valence numbers of protons and neutrons. The latter 
is inspired in the shell model, including explicitly the di- 
agonal two- and three-body residual interactions between 
valence particles and holes. 

In principle the fluctuating part SM also depends on 
the details of the interaction. However, according to 
Strutinsky's energy theorem, the leading contribution 
can be evaluated within the mean field approximation 
which assumes the nucleus is composed of free nucleons 
confined by a one-body potential. It has been shown [ll| 
that even a simple one-body potential, in which the nu- 
cleons are confined inside a spherical rigid sphere (spher- 
ical model from now on), with radius R = rQN nuc x ^ 
(rrj ~ 1.1 fm and N nuc the number of nucleons) describes 
qualitatively some aspects of the experimental SM . How- 
ever, for a more quantitative comparison one has to in- 
clude small multipolar deformations of the spherical cav- 
ity and an effective number of nucleons ^], ^| . While 
the idea of employing a spherical well to describe the in- 
dependent particle model of the nucleus is rather old 0] , 
the corresponding magic numbers, associated with the 
zeros of the spherical Bessel functions, are in only rough 
agreement with the observed nuclear shell closures, even 
when an effective rescaling is employed [ll|. The spheri- 
cal model has shown its best predictive power in systems 
with just one kind of particles, like electrons in spherical 
metal clusters, where shell closures are predicted in close 
agreement with the experimental observation |14|. 

On the theoretical side, a clear advantage of the spher- 
ical model is that SM can be evaluated analytically in the 
semiclassical limit by expressing the exact spectral den- 
sity of a quantum particle in a sphere as a trace formula 
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FIG. 1: (Color online) Mass differences plotted as function of 
Z, A, N, and as an ordered list 



[la | , namely, as a sum over periodic orbits of the classical 
counterpart. In this way explicit expressions for SM are 
available for a nucleus composed of an arbitrary number 
of nucleons. In this letter we will show that this simple 
spherical model is specially suitable for the description 
of the autocorrelations C(q) of SM as a function of the 
number of particles. We will show that a quantitative 
agreement with the experimental mass autocorrelations 
can be obtained without any of the extensions (deforma- 
tions of the sphere and an effective number of nucleons) 
needed for the case of the microscopic contributions to 
the nuclear mass. This is indeed remarkable given the 
simplicity of the model and the complex behavior of the 
nuclear many body problem. 



II. AUTOCORRELATIONS 

Our object of study is the autocorrelation, 
F(q) N 



C(q) 



F(0) N -q 



(2) 



F(q) = J2 SM (i) SM{i + q) 



(3) 



where the sum runs, depending on the case, over the to- 
tal number of nucleons, the neutron number N, or over a 
set includingall possible nuclei as given by the boustro- 
phedon list [HQ]- We shall also investigate C(q) inside 
an isotopic chain, namely, we fix the number of protons 
Z and examine the autocorrelations among all isotopes. 

Autocorrelations are a useful tool in identifying rela- 
tionships between elements in a list or an array. The 
autocorrelation of a constant distribution is also a con- 
stant distribution, and that of a pure harmonic sine or 
cosine distribution will also be an oscillatory distribution. 
On the other hand, the autocorrelation of a random dis- 
tribution is a delta function, peaked at q = plus a small 
random signal for any other q, signaling a null correlation 
length among the elements of the distribution. 

The fluctuating part of the nuclear mass distribution 
can be defined by 



SM, 



exp 



M, 



exp 



M d 



rop i 



(4) 



where M exp is the experimental value for a certain nu- 
cleus according to the nuclear data chart Ij] and M^ rop 
is the prediction of the liquid drop model J5|. As shown 
in Fig 2, the autocorrelation C(q) has a well defined 
oscillatory behavior with clear maxima and minima re- 
lated with the presence of shell closures, as seen in Fig. 
1. When the oscillation amplitude decreases, the posi- 
tion of the first zero in C(q) provides an estimate of the 
size of the region in the nuclear chart where the micro- 
scopic, fluctuating contributions to the nuclear masses 
are strongly correlated. It will be shown that this region 
can include as many of 10 to 15 isotopes or isotones, with 
at least 200 neighboring nuclei significantly correlated. 
The oscillatory behavior of SM is closely related with 
the oscillations in C(q). In what follows it will be shown 
that not only the oscillation length, but also other details 
of these oscillations are well described by the spherical 
model. 

Theoretically SM is expressed as a function of the spec- 
tral density g(E) = J2 t S ( E ~ E i) = S( E ) + S 9( E ) of thc 
one body Hamiltonian (in our case a free fermion con- 
fined in a spherical cavity) as, 



SM = M -M 



with 



N 



M = 2 E i = 2 / E g(E) dE 



(5) 
(6) 



and 



M = 2 / E g{E) dE, 



3 



where E{ are the eigenvalues of the one-body Hamilto- 
nian and g and Sg are the mean and fluctuating part of 
the spectral density respectively. 

The exact (Ep) and smooth (Ep) Fermi energies are 
obtained explicitly as a function of the number of parti- 
cles by inversion of the following relation, 



N„ 



g(E)dE 



g(E)dE. 



(7) 



where N nuc is the number of nucleons (neutrons or pro- 
tons) and the factor two accounts for the spin degeneracy. 
The final expression of SM in term of the spectral density 
is given by, 

SM(N nuc ) = M(N nuc ) - M nuc 

/Ep pEp 
Eg{E)dE-2 I Eg(E)dE. 

In order to compute analytically the autocorrelation 
C(q) we will first evaluate 8M by using the semiclassical 
expression for the fluctuating part of the spectral density 
Sg(E) in a spherical cavity. We then mention how to 
get Ep as a function of the number of particles N nuc . It 
is well known [lfj that, for generic cavities, the smooth 
part of the spectral density g(E) in three dimensions is 
given by, 



9(E) 



2n 2 h 



Vk + S + 



1 

6k 



(8) 



where E = h 2 k 2 /2m, V is the volume of the cavity, S is 
the surface, R\ , i?2 are the radii of curvature and k is the 
scalar curvature. 

For a spherical cavity Eq. JSJl reduces to, 



(E) = ^-E^ 2 R 3 - -R 2 + ^-E- 1 / 2 



3tt 



6ir 



(9) 



In this way the mean Fermi energy is explicitly obtained 
as a function of N nuc by performing the integral in Eq. 
(|7| and then expressing Ep as a function of N nuc . 

In the following section we give a brief account of how 
to evaluate 8g(E) semiclassically by a trace formula in- 
volving only classical quantities. 



III. SEMICLASSICAL EVALUATION OF THE 
SPECTRAL DENSITY IN A SPHERICAL 
CAVITY 

The oscillatory part of the spectral density describes 
the fine structure of the spectrum. These oscillations are 
related with classical periodic orbits inside the cavity 01 
(for an introduction see [Tllll^ ). 

5g(E) = ME) exp(iS a (E)/h + v a ) , (10) 



where the index a labels the periodic orbits, S a is the 
classical action and v a is the Maslov index. As a general 



rule, the amplitude A a {E, L) is a decreasing function of 
the cavity size L but depends strongly on its shape. It 
increases with the degree of symmetry of the cavity. It 
is maximal in spherical cavities and minimal in cavities 
with no symmetry axis. The difference (for the same vol- 
ume) between these two limits can be of orders of mag- 
nitude. 



A. The spherical cavity 

The oscillating part of the spectral density of a particle 
in a spherical cavity of radius R has already been ana- 
lyzed in the literature ITpL IT7j | . Below we provide a brief 
overview and refer to [19] for an account of the details of 
the calculation. 

For a spherical geometry the closed stationary trajec- 
tories are given by planar regular polygons along a plane 
containing the diameter. The length L of the trajectories 
is given by the simple relation L — 2p_Rsin(0) where p 
is the number of vertexes of the polygon and <j> — nt/p 
with t being the number of turns around the origin of a 
specific periodic orbit. Two cases must be distinguished: 
Orbits with p = 2t corresponding with a single diameter 
repeated t times contribute to the density of states as, 



Sg D (E) 



1 



2ttE ^ t 



^2jsm(Aty/E/E ), 



(11) 



where E = 4^ and E = 2 ^' R i ■ For the case p > 2t 
corresponding to regular polygons the contribution to the 
spectral density is given by, 



1 



En V Ei 



££(-1)^(20) 

t=l p>2t 



pir v. 4 



The complete expression for the fluctuating part of the 
spectral density is, 



5g(E) = Sgp(E) + 6g D (E), 



(12) 



where the first term yields the leading correction for suf- 
ficiently large cavities. 

A similar calculation can be in principle carried out 
for a chaotic cavity. In this case the spectral density can 
also be written in terms of classical periodic orbits by 
using the Gutzwiller trace formula. Although an explicit 
expression for the length of the periodic orbits, equiva- 
lent to Eq. (|12fl . is not in general available in this case it 
is still possible to estimate the amplitude of the oscillat- 
ing part by using symmetry arguments. This amplitude 
increases with the symmetry of the cavity. In cavities 
with one or several symmetry axis periodic orbits are de- 
generate, namely, there exist different periodic orbits of 
the same length related by symmetry transformations. 
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It can be shown that the amplitude, as a function of k, 
is enhanced by a factor (kR) 1 / 2 17] for each symmetry 
axis. A spherical cavity has three symmetry axis so the 
symmetry factor S is proportional to S ~ {kR) z / 2 1. 
The factor R is a typical length of the cavity. 

By contrast, chaotic cavities of the same volume have 
no additional symmetries and the symmetry factor S is 
unity, corresponding to the contribution of a single un- 
stable periodic orbit. Consequently finite size effects are 
much more important in cavities with high symmetry. 
We have now all the ingredients to compute the autocor- 
relation C(q) as a function of the number of nucleons in 
the rigid spherical approximation for the nucleus. 

IV. MASS AUTOCORRELATIONS IN THE 
NUCLEAR SPHERICAL MODEL. RESULTS AND 
COMPARISON WITH EXPERIMENT 

In this section we adapt our previous results to the 
specific case of the nucleus. Our aim is to evaluate the 
autocorrelation function C{q) given in Eq. @. We now 
describe the smooth part of the ground state energy M 
by means of the liquid drop model. The fluctuating part 
SM is computed by assuming that nucleons, protons and 
neutrons, are confined in a spherical cavity. Obviously 
this is a mean field approximation that should become 
better as the number of nucleons grows. For comparison 
with the experimental results we will typically remove 
those nuclei with N < 30, a region where the mean field 
approximation is not appropriate. 

In our calculations, the radius R is related to the num- 
ber of nucleons N nuc by R = roNnuc 1 ^ 3 with ro ~ 1.1/m. 
We remark that since neutrons and protons are distin- 
guishable one has to consider these contributions sepa- 
rately, each one with its own Fermi energy but with the 
same radius. We are now ready to write down an explicit 
analytical expression for SM, 

SM = 2 Eg(E)dE-2 Eg(E)dE, (13) 

where the spectral density 5g(E) is given by Eq. (|T21l . 
and Ep is expressed as a function of the number of par- 
ticles N nuc by solving exactly the third order equation in 
Ep, 

with e = Ep/E . Finally the exact Fermi Energy Ep is 
computed by inverting numerically Eq. Q. In all cases 
we assume a mass m p ~ m n ~ 9A0MeV. The sum over 
periodic orbits in Eq. (|12f> has a natural cutoff for scales 
(length of periodic orbits) such that inelastic processes 
which break the quantum coherence are relevant. In or- 
der to account for this fact we have included in the spec- 
tral density Eq. (|T2|l a damping factor k(l) = sill ^f -n* 



where / is the length of the periodic orbit and £ a co- 
herence length that acts as a effective cutoff for I ^S> £. 
Following the estimation of Ref.^J for the nuclear case 
we have set £ ~ 5R. We have checked that the gross fea- 
tures of C (q) do not depend on the cutoff, provided that 
enough periodic orbits are taken into account but other 
parameters like the amplitude of the oscillations of C(q) 
may depend on it. This value of the coherence length can 
be associated with an effective temperature close to 
1 MeV, typical of pairing energies not included in the 
model. 



A. Comparison with experimental results: C(q) as 
a function of the number of neutrons 

We now compute C(q), defined in Eq. J2J), for a nu- 
cleus composed of TV neutrons and Z protons with the 
fluctuating part of the mass given by Eq. lfT3|l . 

First we examine the autocorrelation function as a 
function of the total number of neutrons N. We remark 
that predictions of our model for C (q) are essentially pa- 
rameter free. Since there are many different nuclei with 
the same number of neutrons a proper averaging method 
is needed. In order to proceed C{q) is evaluated as fol- 
lows (see Fig |21 right): we first obtain the analytical pre- 
diction for SM = 6M(N) + SM(Z) for each of the 2140 
combinations of N and Z, then perform an average over 
different nuclei with the same N and finally compute the 
autocorrelation function C(q). The experimental C(q) 
is obtained by using the same averaging procedure. As 
shown in Fig. |3 despite the simplicity of the model, the 
agreement with the experimental results is quite satisfac- 
tory. It accurately reproduces both the amplitude of the 
oscillations and the position of the maxima and minima. 
The agreement between theory and experiment gets bet- 
ter if only heavier nuclei are considered. This is expected 
due the mean field nature of the model. The agreement 
between theory and experiment could be improved if, as 
discussed in [111] , multipolar corrections are considered. 
However we prefer to stick to our parameter-free model in 
order to emphasize that the main features of the autocor- 
relation function are related to the spherical symmetry 
of the problem. We remark that similar results are ob- 
tained if, instead of taking into account all the possible 
combinations of N and Z, we make the simple assump- 
tion SM ~ 25M(N), with r = r {2N) 1 / 3 . 

For the sake of completeness we have also computed 
C(q) as a function of the total number of nucleons A = 
N + Z. As was expected (see Fig EJ a similar degree of 
agreement has been found. 



B. Comparison with experimental results: C(q) as 
a function of the boustrophedon ordering scheme 

By performing averages (for A or N fixed) over the 
nuclear data-chart we may be loosing valuable informa- 
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FIG. 2: (Color online) The autocorrelation C(q) as a function of N = 30 . . . , 154. (Left) and N = 60 . . . , 154 (Right). In both 
cases the agreement with the experimental results (diamond) is quite good. For TV > 60 (Right) it reproduces correctly both 
the amplitude of the oscillations and the positions of the maxima and minima of the experimental data. 
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FIG. 3: (Color online) The autocorrelation C{q) as a function of A = N + Z = 50 . . . , 254. (Left) and A = 116 . . . , 254 (Right). 
In both cases the agreement with experimental results (diamond) is quite good. For A > 116 (Right) reproduces correctly both 
the amplitude of the oscillations and the positions of the maxima and minima of the experimental data. 



tion about nuclear mass correlations. Moreover, since 
cuts along fixed N or A have a small number of nuclei, it 
is difficult to extract definite conclusions. To overcome 
these difficulties, we organize all nuclei with measured 
mass by ordering them in a boustrophedon, namely, a 
ID list composed by 2140 entries numbered as follows: 
Even- A nuclei are ordered by increasing N — Z, while 
odd-A ones follow a decreasing value of N — Z. We have 
evaluated both the experimental and the analytical au- 
tocorrelation C(q), Eq. J2J), as a function of the order 
number of the boustrophedon. For each i = 1, . . . , 2140, 
SM is evaluated by a specific N and Z combination cho- 
sen according to the above classification scheme, as we 



did previously, but in this case we have not performed 
any average. As is shown in Fig0|the agreement between 
theory and experiment is also quite satisfactory for this 
more general correlation function. Both the global os- 
cillatory behavior and the more microscopic details (see 
right plot in Fig. |3J are well reproduced. 

From the above extensive analysis we conclude that 
the main features of the nuclear mass correlations are 
captured by the simple spherical model. As was men- 
tioned previously, our analytical results could be further 
improved by considering small multipole corrections to 
the spherical shape [20|. However it is remarkable that 
our simple spherical model can reproduce in great detail 
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FIG. 4: (Color online) The autocorrelation C(q) as a function of the order number according to the boustrophedon list. (Left) 
Both the result for the spherical model and the experimental data are obtained by considering the 2140 possible combinations 
of N and Z. The agreement between theory (solid line) and experimental data (diamonds) is quite good. (Right) The same 
but now C(q) is plotted only in the window q < 150. 



average properties of the nuclear autocorrelations. 



V. POWER SPECTRUM AND INTEGRABLE 
DYNAMICS 

Finally as a further check of the validity of our results 
we compare the power-spectrum associated to the nuclear 
mass fluctuation SM(i) (i = 1, . . . , N = 2140 is the label 
of the nuclei according to the boustrophedon ordering) 
with the prediction of the spherical model. The discrete 
Fourier transforms of the mass fluctuation is just, 



m 1 E ^l exp « 

~fN ^ <TrmB \ N J 



(15) 



with the root-mean-square (rms) deviations given by 

1/2 



1 N 



(16) 



where 5M(j) is either the experimental or the analyti- 
cal fluctuating part of the nuclear mass. The decay of 
the associated power spectrum S(k) = |F(fc)| 2 provides 
information about the type of dynamics of the model. 
Thus it can be shown [21j that, for scales roughly in 
between the shortest periodic orbit and the mean level 
spacing, a power law decay S(k) ~ k~ a with a = 4 cor- 
responds to integrable classical dynamics. In Fig [5] we 
observe a close agreement between the power spectrum 
of the spherical model and that of the experimental nu- 
clear masses. Moreover the decay in the range ~ [1,3], 
which includes frequencies between those associated with 
the mean level spacing and with the shortest periodic or- 
bit, follows a power-law with a ~ 4.2, in agreement with 
the prediction for integrable dynamics. Based on these 



results we suggest that the power spectrum could be uti- 
lized as an effective test to check whether a strongly inter- 
acting many body system is indeed close to integrability 
or not. 

The power spectrum of differences between measured 
masses and those calculated in different models has been 
studied in Q- A gradual vanishing of the slope a was 
observed as more sophisticated and realistic models were 
utilized. For the most realistic models a white noise a — 
(all frequencies have equal weight) signal was found. 
For a detailed study of intermediate situations we refer 

to m. 
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FIG. 5: (Color online) Power spectrum S(k) of the nuclear 
mass fluctuation for the boustrophedon ordering. As is ob- 
served, the agreement between theory and experiment is very 
good in the low and intermediate frequency region. Their 
power-law decay with a ~ 4.2 is close to the result a = 4 
predicted for classically integrable systems. 
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VI. CONCLUSIONS 

A simple semiclassical analysis, where protons and 
neutrons are described by free particles bouncing elas- 
tically, back and forth inside a rigid sphere, has been 
shown to nearly reproduce the autocorrelations of the 
differences between measured nuclear masses and those 
calculated using the liquid drop model. The results are 
remarkable, offering a different insight on the microscopic 
corrections needed to describe nuclear masses with pre- 
cision, ft also has been shown that it is possible to per- 
form autocorrelation analysis of nuclear mass differences 



along very long chains of isotones, isotopes, isobars and 
other chains, a task generally considered very difficult to 
perform [2]|. While interesting in themselves, these re- 
sults could also provide a theoretical explanation of the 
amazing success of the two dimensional Fourier analysis, 
performed in the Z — N sp ace, in the description and 
prediction of nuclear masses |24J . 
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